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ON SUBSONIC FLOW PAST A PARABOLOID OF REVOLUTION 
By Carl Kaplan 


SUMMARY 


The Janzen-Rayleigh method of* expansion in powers of the stream 
Mach number M^ is utilized for the calculation of the velocity poten- 
tial for steady subsonic flow past a paraboloid of revolution. Only the 
first two terms of this expansion are calculated, the first term being 
the incompressible expression for the velocity potential and the second 
being the term in M^ 2 . A closed expression is obtained for the second 
term in the form of a double infinite integral which contains Bessel 
functions under the integral signs. The methods of evaluating such inte- 
grals are not very numerous. Unfortunately, the present integral does 
not yield to any of them. No attempt is made in the present paper to 
evaluate numerically this double infinite integral. Expressions for 
the fluid velocity axe given in the form of correction factors by which 
the corresponding expressions for incompressible flow are multiplied in 
order to take into account the effect of compressibility. 


INTRODUCTION 


The Janzen-Rayleigh method for the calculation of subsonic flow 
past an obstacle, by expansion of the velocity potential in powers of 
the undisturbed stream Mach number, has been extensively applied to two- 
dimensional problems. A considerable void, however, exists in the liter- 
ature insofar as applications of this method to three-dimensional axisym- 
metric problems are concerned. Except for the case of flow past a sphere 
and sporadic efforts to treat the next simplest case of axisymmetric flow 
past a prolate spheroid, little has been accomplished in this area of 
subsonic compressible- flow theory. Clearly, then, future work utilizing 
the Janzen-Rayleigh method can be expended profitably on axi symmetric- 
flow problems. A first step in this direction has been taken by 
A. L. Longhorn who recently, at the suggestion of M. J. Lighthill, recon- 
sidered the problem of subsonic flow past a prolate spheroid. (See 
ref. 1.) The problem treated in the present paper was begun sometime 
before the author became aware of the results of Longhorn. Fortunately, 
the two problems complement one another in the sense that the one treated 
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by Longhorn is limited to closed bluff bodies, whereas the present one 
is concerned with a family of semi-infinite elongated bodies with round 
noses . 

The choice of the paraboloid of revolution as the solid body was 
made for several reasons. First, the paraboloid being a semi-infinite 
body and, in fact, a limiting case of the prolate spheroid, it was 
thought probable that the analysis might involve functions of a more 
elementary nature than for the case of a closed body. Second, the fluid 
speed at the surface rises monotonically from zero at the stagnation 
point to the undisturbed-stream value at infinity. Therefore, the criti- 
cal value of the stream Mach number is unity, and hence there can be no 
transonic influence in the entire subsonic range. Finally, the Janzen- 
Rayleigh method being a thoroughly reliable one, the present investiga- 
tion should provide useful information with regard to the question of the 
accuracy of the small-disturbance method for the calculation of compres- 
sible flow past slender bodies in the neighborhood of the stagnation 
point . 


The Janzen-Rayleigh method has often been described in the liter- 
ature. In the present paper, therefore, only those equations necessary 
for the formulation of the problem are used. The problem to be considered 
is the subsonic flow past a paraboloid of revolution fixed in a uniform 
stream of velocity U in the negative direction of the axis of symmetry. 
The nature of axisymmetrical flow is such that the motion is the same in 
every (meridian) plane through the axis of symmetry. The position of a 
point in a meridian plane may be fixed by rectangular Cartesian coordi- 
nates x,y with the origin at the focus of the parabolic meridian pro- 
file. (See fig. 1.) With the radius of curvature at the nose as the 
unit of length, the equation of the meridian profile becomes 


For steady subsonic axisymmetric flows, the equation satisfied by 
the velocity potential <f is the equation of continuity 
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( 1 ) 




where 
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(3) 


NACA TN 3700 


3 


and 


q nondimens ional speed of fluid with U as unit of velocity 

p density of fluid moving with speed q 

density of fluid moving with undisturbed speed U 
Mach number of undisturbed stream, — 

C oo 


c^ speed of sound in undisturbed stream 

7 ratio of specific heats at constant pressure and constant 

volume 

If, now, a series expansion for 0 in powers of is assumed, 

then 


n =0 


0X 


2 n 


where the first term 0 Q is the velocity potential in incompressible 
flow. In this paper, only 0^ is calculated so that the form of 0 
to be found is 


0 = 0 O + Mco 2 0! (V 

Thus, this expression for 0 is substituted into equation (2), and equa- 
tion ( 3 ) is utilized with 



2 

Then, when the coefficients of the zeroeth and first powers of are 

equated to zero, the following pair of equations is obtained: 



J (6) 


where q 0 is the local velocity in incompressible flow. 
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The velocity potential for incompressible flow past the paraboloid 
of revolution defined by equation (l) is readily found to be 


0 o = -x + i log^x + VJ77) (7) 

where 0 O is nondimensional with the undisturbed velocity U as unit 

of velocity and the radius of curvature at the nose as unit of length. 
Equation ( 7 ) satisfies equation (5) and the boundary conditions of van- 
ishing normal velocity at the surface and of vanishing disturbance veloc- 
ity infinitely far from the paraboloid. The expression for the magnitude 
of the incompressible fluid velocity is 


<o 2 -l 




1 - 


x 


+ y c 


2y* 




(8) 


+ y* 


At this point it is convenient to introduce a new set of independent 
variables. Ideally, the appropriate coordinate system for the present 
problem is a parabolic one which defines mutually orthogonal families of 
confocal parabolas in a meridian plane. Thus, the conformal transformation 


= x + iy = (| + iT )) 2 = 


gives 


x = 




2 


y = 2 |tj 

J 

The elimination of the variable T) in equations ( 9 ) yields 


(9) 


X - 



Then the surfaces 4 = Constant = | Q are confocal paraboloids of revolu- 
tion with the focus at the origin and the radius of curvature at the nose 
equal to 2| Q 2 . With this radius of curvature as unit of length, equa- 
tions ( 9 ) are nondimensional and the solid boundary is given by ^o = ^ 

or equation (l). 
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The equations corresponding to equations (5) and (6) are 

U* W) + w) - 0 (10) 


and 


^i\ , b_L = 1 i_ 

8a\ 8a J d(3 \ P ap / 2 da 


“ - # 


1 8_ 

2 8p 


P (-Jo 2 - l ) 


80 o 

8T 


( 11 ) 


where 


|2 = a and t]^ = p 


Also, equations ( 7 ) and (8) become, respectively, 

(ZS o = -a, + P + i log 2a 


and 




a + P 4a(a + P) 

Then, equation (ll) for the velocity potential 0-^ takes the form, 

*-(.£)♦» (3) 


(12) 


(13) 


8a\ 8a / 8pV 8p / 2 ( a + p) 2 2 ( a + P) 2 l6a 2 (a + P) l6a(a + p) 2 

(14) 

The complementary solution of this partial-differential equation is 
easily found by assuming a solution in product form: 


= A(a)B(p) 


Then the homogeneous form of equation (l4) becomes 
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or 


da \ da/ 


_d_ 

dp 


( e 




(15) 


where A is an arbitrary real constant. Equations (15) are equivalent 
forms of Bessel's differential equation for cylinder functions of order 
zero. Therefore, 


A = C 0 (iA\/a) 


and 


B = C 0 (Up) 


where C Q represents any cylinder function of order zero. 

The nature of the present problem is such that the cylinder function 
involving the variable a shall behave well for a >°° but not neces- 
sarily for a — > 0 since a = ^ at the solid boundary and never takes 

on values less than that. On the other hand, the cylinder function, 
involving (3 must be of such type that it behaves well for the entire 
range 0 ^ (3 ^ 00 . Clearly then. 


A = K 0 (A^a) 


and 

B = J 0 (Afp) 

where K 0 is the modified Bessel function of the second kind and zero 
order and J Q is the Bessel function of the first kind and zero order. 
The graph of K 0 resembles a rectangular hyperbola in the first quadrant 
and that of J 0 , a damped cosine wave. The general complementary solu- 
tion of equation (l4) can then be written as 

OO 

- y* c n^o (\i^)^o(\i^) 
n=l 

where the values of c n are arbitrary constants to be determined by 
means of the boundary condition that the fluid-velocity component normal 
to the surface of the paraboloid shall vanish. 
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A suitable particular integral of equation (l4) is now sought. 

Unf ortunately, the right-hand side of that equation cannot be separated into 
a sum of products of functions of a alone and of P alone; hence, the 
task of finding a particular integral directly from equation (l4) is prac- 
tically futile. However, by temporary use of polar coordinates as inde- 
pendent variables the desired separation of variables can be achieved. 

Thus, 


x = r cos 6 
y = r sin 0 


(17) 


where the radius vector r is nondimens ional with the radius of curvature 
at the nose of the paraboloid as unit of length and where the angle 0 is 
measured positive counterclockwise. (See fig. 1.) 


The polar equation of the parabolic meridian profile is 


1 + cos 0 


( 18 ) 


and the expressions for the incompressible velocity potential and fluid 
speed given by equations (7) and (8), respectively, are 


0 O = -r cos 0 + i log r(l + cos 0) 

a 2 = p _ 1 + 1 

r 2r^(l + cos 0) 


(19) 


Thus, at the upper surface of the paraboloid of revolution, the expression 
for the incompressible fluid speed is simply 


q = sin i 0 
^o 2 

It is interesting to note the curious fact that this expression for the 
velocity at the surface of a paraboloid of revolution is precisely the 
same as that for the velocity at the surface of a two-dimensional para- 
bolic cylinder. The explanation lies in Munk's rule which states that 
the surface velocity on any ellipsoid immersed in a uniform flow along 
a principal axis is the projection of the maximum velocity (in the pres- 
ent case, unity) on the tangent plane to the surface. (See ref. 2.) This 
is obviously the same for both a parabolic cylinder and a paraboloid of 
revolution in view of the fact that Munk's rule includes the two- 
dimensional case. 
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When equations ( 17 ) and ( 19 ) are utilized, differential equation ( 6 ) 
for the velocity potential 0 1 takes the following form: 


( 2 

. d_ 

"(1 u2)^ 1_ 

__ln , 1_ _ 

I 

OJ 

H 

i — 1 

1 

V dr / 

dp 

\-L - M' 

d[i 

2 2r 8r 2 

1 + p (i + p) 2 


8_ 

dr 


where \i = cos 0 • 


( 20 ) 


Note that the desired separation of variables has been achieved on 
the right-hand side of equation (20). The task of finding a particular 
integral then becomes a routine problem. The following particular inte- 
gral 0 1 of equation ( 20 ) has been constructed with the point in mind 

that it be well behaved everywhere at the surface of the boundary and in 
the field of flow: 


0 lp ’ ^ & 108 r(l + 11)5 + ^2 


1 - log 2 - 


1 + p 


— p log 2 (l + p) + ^1 + ^ p log 2 ^ log(l + p) - 


i r ^ iog( l + p) - log 2 

I p 

d 1 0 1 - M- 


( 21 ) 


This particular integral may be expressed in terms of parabolic coordi- 
nates by means of the following equations of transformation: 


r = a + (3 

a - P 


a + P 


The general solution for the velocity potential 0^ then becomes 




1 a - p 1 lQg 8a3 


1 k a + p 8 a + p ( a + p ) 2 8 (a + p ) 2 

fl + - - ~ log 2) log - 
\ 2 a + P / a 


— - log 2 - - 

2 2 a 


la- P . 2 2a 

— - log + 

4 a + P a + P 

a-p 

1 a - p r^Tp log(l + h) - log 2 ^ 

2 a + PJq 1 “ M- 


2 a 


+ P 


+ 2_ c 

n=l 


n K o(vk) J °M) 


( 22 ) 


( 23 ) 
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where the arbitrary constants c n are determined by means of the bound- 
ary condition that the normal component of the fluid velocity at the sur- 
face of the paraboloid shall vanish. In addition, there is the require- 
ment that the disturbance velocity shall vanish at infinity for points 
not near the paraboloid. Now, 0 O , being the solution for incompressible 

flow, satisfies these boundary conditions, and hence 0^ must also sepa- 
rately satisfy them. An examination of equation (23) shows that the dis- 
turbance velocity vanishes at infinity. The normal boundary condition 

/d0,\ 

— - = 0 then yields the following equation for the determination of 

\da /(x=7j. 

the arbitrary constants: 


11 Vo(\#) = 

n=l 


(1 + 2 p ) 5 

(1 + 2p)log 2 

I log 2 (l + 23) + 1(3) 


|(1 - 23) [-1 + log( 1 + 23)] + 
l - 43 ~ 


(l + 23) 4 


- | log 2 2 + 


5 


= f(3) 


where , for convenience of expression. 


d 


n 



(24) 


and 


i(p) 


1-23 

1+23 


log( 1 + p) - log 2 


0 


1 - p 


dp 


Now, the range of the variable 3 along the upper surface of the solid 
boundary extends from 0 to ». Therefore, the determination of the 
complementary function 0^ c must proceed along lines corresponding to 

the passage from a Fourier series to a Fourier integral. Thus, multi- 
plying both sides of equation (24) by /3Jo(\i\/P) (a ^P and integrating from 
0 to an arbitrary upper limit b give 
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r b 

dn = ^- 7 ^ F(p)\fal 0 (A n #)d# (25) 

[bJi(bA n )J J ° 

where A^_, T\2> • * * are different values of A for which J 0 (bA n ) = 0. 

In equation (25), use has been made of the well-known orthogonality- 
condition 



(\a ^ ^n) 


(\i - 7v n ) 


The quantity is now defined as bA n /b, and the upper limit b 

is taken to be very large. The difference between two consecutive values 
of 0^ may then be obtained by considering large zeros (bA n ) of J Q * 
Thus, the approximating formula for the nth zero is 


(26) 

Hence, the distance between consecutive zeros approaches rt and the 
difference between two consecutive values of becomes 



or 



1 

b 



Now, for large values of bA n , 

sin 

Jp(bA n ) = 

or with the aid of equation (26) , 


i n (bA n ~ f) 

\j| bA n 


Jl(bAn) 



MCA TN 3700 


11 


It follows then from equation (25) that 


*n 


cnbcn / F(p)\/pJ 0 (cjDVp)d.7p 

Jo 


Introducing this expression for d n into equation (l6) for the comple- 
mentary function 0q c > replacing the summation by integration, and letting 
the upper limits go to infinity lead to the following expression: 


Jo Mi 

Inserting this expression for the complementary function into equation (23) 
then yields the exact form for the second term in the Janzen-Rayleigh 
method for the calculation of subsonic flow past a paraboloid or revolu- 
tion. Equation (27) can be verified by means of the boundary condition 




= F(P) 


and the recurrence relation 

K (oj'/a ) 


da 


—1 


= --^ w (*) 


v/2 H1/2 


Thus the well-known Fourier-Bessel integral (ref. 3) is obtained: 

poo pOO 

F(p) = / CO dco / F(t 2 )J 0 (cot)J 0 (co\/p)t dt 

Jq Jn 


NUMERICAL CONSIDERATIONS 


Attempts to evaluate the right-hand side of equation (27), or even 
to reduce it to a single infinite integral, have thus far proved fruit- 
less. It appears certain that, in order to obtain numerical results, 
equation (27) for 0-j_ c or its derivatives with regard to a and 3 

must be calculated directly. This calculation is not done in the present 
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paper. Rather, in preparation for such numerical computations it is 
necessary to evaluate the integral appearing in the expression (eq. (24)) 
for F(p); thus, 

1-2(3 

I(P) g r 1+2(3 log(l + fj) - dtx (O^P^co) 

J 0 1 “ M- 


Then, the following expansions are considered: 


“ „ n 
log(l + n) = -51 (-1) Hj- 


n=l 


and 


(-i) n 

log 2 -51 


n=l 


n 


It follows that 


log( 1 + n) - log 2 

1 - p 


/ n ^n n 

(-1) ^ — m-1 


n=l n m=l 


Hence 


w ( t \ n n m 

I(p) = \ ( - 1 I V- 


n=l n m=l m 


(-1 £ p < l) (28) 


Note that the convergence has been improved by the act of integration to 
include p = -1. Thus, 


00 , ,n 

lim l(p) = - S, 


P — * 0 

or 

M- 


n=l 


n n 


where 


s -^1 

n ~bi m 


Now 


log 2 2 . ± Ail! s „ . 


00 , . n 

_ y (- 1 ) 

n=l n n bl n2 
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and 


Therefore, 


n=l n ' 


(-D n __*£ 

2 12 


lim 1(3) = i log d 2 - 2- » - 0.58224 
0— >0 12 

For 0 — >00 or p — >-l, equation (28) gives 

" v n n , _ ,m 


- «»>-z4Fr4> 

0— >°° n=l m=l “ 

or after rearrangement of terms on the right-hand side 


( 29 ) 


lim l(|3) =■ 


(_l) n 

S_ + 


n=l 


n 


n 


n ^=1 2n(2n + 1) 


= - lim 1(0) + XT - 

0— >0 n=l 2n(2n + 1) 

The series on the right-hand side can be rendered more rapidly convergent 
by repeated application of Kummer's transformation (ref. 4); thus, 

lim l(0) = i log 2 2 + — + — 2Z — + 

0 — >oo 12 2 60 4 n=1 (n + l)(2n + $)(2n + 5) 


15 


^n 

^1 n(2n + l)(2n + 3)(2n + 5) 


« 1.060 

Finally, for 0 < 0 < °° or p 2 < 1, 


1 ( 3 ) = 


/ , x n n w 

(-J-) yr 

-1 n m 

n=l m=l 


°° n+1 n 

= log 2 log( 1 - p) - 21 


>=i n+1 ■ 


n=l 


m=l 


m 


00 2n+l 

= log 2 log(l - M-) - 2 l0g ^ 1 + V) 10 g(! - H-) + 51 g - n - t 2n 

n=l 
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where 


t 2n = 1 




S 


n 


Thus, 


1(0) = J log — ^ 

2 1 + 20 


log 2(1 + 20) + 


'2n 


Zn 2n + 1\1 + 


n=l 


1 - 


20\ 

20 / 


2n+l 


(0 < 0 < oo) 


(30) 

With the aid of equations (29) and (30) the function F(0) can be com- 
puted for all finite values of 0. For 0 < 0 < 00 , the expression for 
F( 0) in equation (24) becomes 


F(0) = 


(1 + 20) 5 


|(1 - 20) [-1 + log(l + 20)] + (1 + 20) log 2 


1-40 

(1 + 20) 4 


i log 20 log 2(1 + 20) + y~ 
^ n=l 


t 2n /I - 20\ 

2n + 1\1 + 20/ 


2n+l 


Table I lists values of F(0) and l(0) for the range 0 < 0 < 9, 
and figures 2 and 3 show the graphs of each. For the later purpose of 
calculating the velocity distribution in the neighborhood of the stag- 
nation point and also along the axis of symmetry of the flow, general 
expressions are' presented in the following section. 


CALCULATION OF FLUID VELOCITY AT THE SOLID BOUNDAKY 
AND ALONG THE AXIS OF SYMMETRY OF FLOW 


The fluid velocity at a point on the body is given by 


where 

and 


ds i 


q 

S S J 


is the element of arc in the 


d0 

direction in which 


T] 


increases 


t 2 _ dz dz 

dl 

Then, in terms of the variable 0 = tj2. 
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(3D 


With the aid of equations (12) and (23 ) , equation (31) yields the fol- 
lowing expression for the fluid speed along the solid boundary for the 
range 0 ^ p <1: 


1 sM 00 ^SG| + 


1 - P 


log 


log 


(1 + 2f3)H 1 + 2P 2 (! + 2(3) 


+ 21 ( (3) 


1 - 2|3 log(l + 2p) _ (1 - 2(3) (l + 4(3) lQg 2 + 1 

(1 + 23)^ 4p 2(1 + 23)^ 1+2(3 2(1 + 2P) 2 


6-3 log 2 

2(1 + 2 p ) 5 


1 

(1 + 2 ( 3 )^ 


log 2 + 



(32) 


where q = 1/ — — — 

SO V 1 + 23 

face in incompressible 


is the magnitude 
flow and c 


of the fluid velocity at the sur- 
is given by equation (27) • 


by 


The fluid velocity along the axis of symmetry of the flow is given 


qx " It ~ (aS) p=0 


Then with the aid of equations (12) and (23), 


^xM„ ^xo 


1 + 


KT 

4x(2x - 1) 


1 - log 8x + i(l i- f§) + 

K 7 p=o 


(33) 


where q = -1 + — is the velocity of the fluid along the axis of 
xo 2x 

symmetry in incompressible flow. Equations (32) and (33) have the form 
of correction factors by which q go and q xo are multiplied in order 

to take into account the effect of compressibility. It may be desirable 
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to utilize the arc length s p of the parabolic meridian profile as the 
reference variable rather than 3 . The arc length s^, with the radius 
of curvature at the nose as unit of length, is given by 


or 


= i |\/ 23 ( 1 + 23 ). + log(V 23 + \Jl + 23 


So = 


n 0 

sin n 1 + sin 75 

-| — los ■ — — 

1 + cos 9 2 CO s 


(34a) 


(34b) 


l 


For purposes of comparison, consider the two-dimensional case of uniform 
subsonic flow past a parabolic cylinder. The following formulas, corre- 
sponding to equations ( 32 ) and (33)> 3-re easily obtained from reference 5 



r 



OaML = %o ‘ 

co S 2 | 

2 + cos |Y 9 ? QS e 9 - 4 cos | log 2 cos |j 




\ sin 2 - 

J 


J 

(35) 


and 


^xMqq Sco 


1 - i M 2 

p 00 


€ 


log | X 


where q = sin - and 


^2x 4x(*/x - /2) 4x(\fx - \/2)_ 

1 


(36) 


= -l + 


y/2x 


CONCLUDING REMARKS 


In conclusion, the present paper provides an additional example to 
the sparse literature on subsonic axisymmetrical flow . The attempt has 
been made to choose as solid boundary a shape which does not require 
involved and cumbersome analysis but, at the same time, which is of 
interest to both theoretical and applied aerodynami cists. Attention is 
particularly directed to the double infinite integral of equation ( 27 ). 
Such integrals, involving Bessel functions under the integral signs, are 
not only of great interest to the pure mathematician but are also of 
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extreme importance in many branches of mathematical physics. The addi- 
tion of another one to the large number of such integrals which have 
already been evaluated can generally be counted upon to aid in the solu- 
tion of many problems in varied fields. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 14, 1956. 
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TABLE I 

TABULATED VALUES OF THE FUNCTIONS l(p) AND F(p) 


p 

i(P) 

F(p) 

P 

i(P) 

F(p) 

0 

-0.58224 

-1.62932 

0.875 

0.20199 

0.07033 

.02083 

-.54184 

- 1.18932 

1.000 

.25103 

.04879 

.04167 

-.50378 

-.85251 

1.200 

.31729 

.02452 

.06250 

-.46788 

-.59387 

1.400 

.37215 

.00895 

.08333 

-.43393 

-. 39^90 

1.600 

.41845 

-.00105 

.10417 

-.40176 

-.24175 

1.700 

.43900 

-.00461 

.12500 

- . 36666 

-.12398 

1.800 

.45812 

-.00746 

.14583 

-.34222 

-.03368 

1.900 

.47591 

-.00974 

.16667 

-.31459 

.03524 

2.000 

.49255 

-.01155 

.18750 

-.28825 

.08742 

2.500 

.56184 

-.01612 

.20833 

-.26309 

.12650 

3.000 

.61442 

-.01693 

.22917 

-.23905 

.15527 

3.500 

.65589 

-.01633 

.25000 

-.21602 

.17595 

4.000 

.68955 

-.01523 

.31250 

-.15249 

.20490 

5.000 

.74108 

-.01283 

•37500 

-.09606 

.20510 

6.000 

.77160 

-.01074 

.43750 

-.04555 

.19195 

7.000 

.80804 

-.00904 

.50000 

0 

.17329 

8.000 

.83123 

-.00771 

.56250 

.04132 

.15309 

9.000 

.85016 

-.00664 

.62500 

•75000 

.07901 

.14535 

.13333 

.09829 

00 

1.06000 

0 



Figure 1.- Profile of paraboloid of revolution in meridian xy -plane. 
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Figure 2.- The boundary function F(p) plotted against p. 
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Figure 3 .- The integral I ((3) plotted against 3 . 
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